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Abstract: Given a quadratic two-parameter matrix polynomial Q(X, fx), we de- 
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parameter eigenvalue problem. 

Key words: matrix polynomial, two-parameter matrix polynomial, quadratic 
two-parameter eigenvalue problem, two-parameter eigenvalue problem, lineariza- 
tion 

AMS classification: 65F15, 15A18, 15A69, 15A22 

1 Introduction 

We consider two-parameter quadratic matrix polynomials of the form 

Q(X, fj) = A + XB + fiC + X 2 D + XixE + y?F (1) 

where A, \x are scalars and the coefficient matrices are real or complex matrices 
of order n x n. If (A, fi) G C x C and nonzero x G C n satisfy Q(X, fj,)x = 0, then 
x is said to be an eigenvector of Q(X, /i) corresponding to the eigenvalue (A, /i). 
The classical approach to solving spectral problems for matrix polynomials is 
to first perform a linearization, that is, to transform the given polynomial into 
a linear matrix polynomial, and then work with this linear polynomial (see 
[TB| [HJJ [T71 [2H EI] and the references therein). Therefore, given a quadratic 
two-parameter matrix polynomial Q(X, /j), we seek linear two-parameter matrix 
polynomials L(X, fx) = XA± + fiA 2 + A 3 , called linearizations, which have the 
same eigenvalues as Q(X,fi). 
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The one-parameter matrix polynomials have been an active area of research 
in numerical linear algebra J2TJ 123 [28] . In [21], Mackey et al. have investi- 
gated the one-parameter polynomial eigenvalue problem extensively and they 
have produced vector spaces of linearizations for a one-parameter matrix poly- 
nomial by generalizing the companion forms of the one-parameter polynomial. 
Adopting a similar approach we derive a set of linearizations of a quadratic 
two-parameter matrix polynomial. 

The quadratic two-parameter eigenvalue problem is concerned with finding 
scalars A,/i6 C and non-zero vectors X\ G C ni , x 2 G C™ 2 such that 

f Q 1 (X,fi)x 1 = (Ai + XB 1 + fid + X 2 D 1 + XfiEt + fx 2 F 1 )x 1 = . , 

\ Q 2 (A, /j)x 2 = (A 2 + XB 2 + aC 2 + X 2 D 2 + XuE 2 + fi 2 F 2 )x 2 = [Z) 

where Aj, B i: . . . , Fi G C niXni ,« = 1,2. A pair (A,/i) satisfying (j2]) is called 
an eigenvalue of (T5]) and Xi <E> x 2 , where ® is the Kronecker product, is the 
corresponding eigenvector. This problem appears in stability analysis of dif- 
ferent systems, for example, time-delay systems of single delay [HI [TOj [131 121] • 
The standard approach to solving ([2]) is by linearizing the problem into a two- 
parameter eigenvalue problem of larger size and then by converting it into an 
equivalent coupled generalized eigenvalue problem which is then solved by nu- 
merical methods, see [2U E] . 

Given (T5|), we seek a two-parameter eigenvalue problem 

/ L 1 (A,/i)w 1 :=(A( 1 )) + A J B( 1 )+/iC( 1 ))w 1 = , , 

\ L 2 (A, fj,)w 2 := (A( 2 )) + XB^ + iuCW)w 2 = [6) 

with the same eigenvalues, where Wi G C 3 "^ \ {0} and AW, B®, G C 3niX3rii , 
i = 1, 2. In such case ([3]) is called a linearization of (J2|). 

The choice of a linearization may have an adverse effect on the sensitivity 
of the eigenvalues. Therefore, it is important to identify potential lineariza- 
tions and describe their constructions. In this paper, we develop a systematic 
approach that enables us to generate a class of linearizations for a quadratic 
two-parameter eigenvalue problem. 



2 Linearizations for quadratic two-parameter 
matrix polynomial 

In this section we construct a set of linearizations of a quadratic two-parameter 
matrix polynomial. 

Definition 2.1 A In x In linear matrix polynomial L(A,/i) = AAi + 

fiA 2 + A 3 is a linearization of an n x n matrix polynomial Q(X, a) if there 
exist polynomials P(X, /i) and R(X, /i), whose determinant is a non-zero constant 
independent of X and fx, such that 







l (l-l)n 



P(\,h)L(\,vl)R(\,(jl). 
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Let Q(X, //) be a quadratic two-parameter matrix polynomial given by 



Q(\, p) = X 2 A 20 + n 2 A 02 + \fiA n + XA W + fxA 01 + A 



oo 



where the coefficient matrices are of order nxn. Assume that x is the eigenvector 
corresponding to an eigenvalue (A,/i) of Q(X,fx), that is, Q(X,fJ,)x = 0. With 
a view to constructing linearizations of Q(X,/i), we denote x = xoo,Ao;oo = 
x w ,fix 00 = x i- Then we have 

A 20 (Xx 10 ) + A 02 (fix 01 ) + A n (Xxoi) + A w x w + A 01 x 01 + A 0Q x 00 = 0. (4) 

Consequently we have 
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the eigenvector corresponding to an eigenvalue (A,//) of Q(A,/x) if and only if 
L(X, n)w = where w = A g) x and L(X, /i) = XAi + /M 2 + A 3 , that is, w is 
the eigenvector corresponding to an eigenvalue (A,/i) of L(X,/j>). We show that 
L(A,/i) is a linearization of Q(X,[i). 
Define 



£(A,m) 

F(A,/i) 



XI n 
/J,I n 
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In 









In 





„ fiA 02 + XA n 
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01 



XA 20 + A 10 

In 





Notice that E, F are unimodular quadratic two-parameter matrix polynomials, 
that is, determinants of E and F are constants. Then we can easily check that 



F(\,n)L(\n)E(\,n) 



Q(X,fz) 






hn 



(6) 



Thus we have detQ(A, fj,) = 7 det L (A, fi) for some 7 7^ 0. This implies L(X,fi) 
preserves the eigenvalues of Q(X, fx) and hence is a linearization of order 3n x 3n. 
We call this linearization the standard linearization of Q(X,fi). It is interesting 
to observe that the linearization proposed in [24] is up to some permutations of 
block rows and columns of the standard linearization. 
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However, for the standard linearization we have 



L(X, n) ■ (A ® x) = [(Q(X, fi)x f ... of for all x G C n , 



(7) 



and therefore, any solution of (jSJ) gives a solution of the original problem 
Q(X, fi)x = 0. Further, by (j7|) we have 





"A/„" 




"Q(A,/x)" 




"l" 


L(A,/x)-(A«)I n ) = Li(A, /x) 









= ei<g>Q(A, /x), ei = 







/„ 













We restrict our attention to the equation (jSJ) which is satisfied by the stan- 
dard linearization. It would be worthy to find linear two-parameter matrix 
polynomials L(A, /x) that satisfy 



L(X, n) ■ (A ® I n ) = v ® Q(X, /x) 
for some vector t> G C 3 . Therefore, we introduce the notation 

{v® Q(A,/x) : u G K 3 } 



(9) 
(10) 



and define 
L(Q(A,/x)) : 



L(A,/x) 



AAi + /xA 2 + M, Ai G K 3 " x3n : L(X, zx) • (A ® J n ) g Vq 

(11) 

Note that L(<5(A,/x)) 7^ as the standard linearization L(A,/x) G L(<5(A,/x)). It 
is easy to check that L(Q(A, /x)) is a vector space. If L(A,/x) G L(Q(A,/x)) for 
some v G C 3 then call v is an ansatz vector associated with L(A,/x). To investi- 
gate the structure of each L(A,/x) G L(Q(A, /x)), we define a "box-addition" for 
three 3n x 3n block matrices as follows. 

Definition 2.2 Let X,Y, Z G C 3nx3n fre i/iree fr/oc£; matrices of the form 



X 
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X12 


X13 




Y u 


Y u 


Yu 




Z n 
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where '+ 7 is the usual matrix addition. 
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For the standard linearization L(A, /x) = XAi + fiA 2 + A 3 e L(Q(A,/i)) we 
have 
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= ei ® [A 20 A u A 02 A 10 A 01 
Thus we have the following lemma. 



4o] . 



Lemma 2.3 Let Q(X, /j) = \ 2 A 20 + ^ 2 A 02 + X/iA u + XA W + /iA 01 + A 00 be 
a quadratic two-parameter matrix polynomial with real or complex coefficient 
matrices of order nxn, and L(X, ji) = XAi + /iA 2 + A 3 a3nx3n two-parameter 
linear matrix polynomial. Then 



L(X, /i) • (A <g) 



v <g) Q(X,/i) 

Ai ffl A 2 ffl A 3 = v <g> [A : 



20 



An 4)2 4o 
Proof: Computational and easy to check. 

Example 2.4 Consider a quadratic two-parameter polynomial 

Q(X, /i) = X 2 A 20 + fi 2 A 02 + X/iAn + XA 10 + fiA 01 + A 



A i A 00 ] . 



00 



where 4j e C nxn and 



A 20 A u + A 20 
L(X, /i) = A A 20 A 00 

2A 20 A 02 + 2A n 

—A 2 o Aq 2 
+/i A u - Aqq A 02 
— A 02 2A 02 

Then L(A,/i) G L(Q(A,/i)) since 

A 1 ffl A 2 ffl A 3 = [1 1 2] T <g> [A 
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Using Lemma [231 we characterize the structure of any L(X, fi) G h(Q(X, [/,)). 



Theorem 2.5 Let Q(\, //) = X 2 A 2Q + [± 2 A m + X/j,A n + XA 10 + fiA 01 + A 00 be 
a quadratic two-parameter matrix polynomial with real or complex coefficient 
matrices of order n x n, and v G C 3 . Then a linear two-parameter matrix 
polynomial L(X, /z) G L(Q(A, fi)) corresponding to the ansatz vector v is of the 
form L(X, fi) = XAi + fiA 2 + A 3 where 



A, 


= [vi 


§» A 20 


-Y,+ 


v®A u -Z x + v® A 10 ] 


A 2 


= [Yi 


v <g 


A 02 -Z 2 + v ® A 01 ] 


A 3 


= fa 


z 2 


v ® A 00 ] 






~Y U ~ 




'Zn 




Z\ 2 




where Y\ = 


Y21 


,Zi = 


Z 2 \ 


,z 2 = 


z 22 


G C 3nxn are arbitrary. 




Y 31 




Z31 




Zz 2 





Proof: Let M. : L(Q(A, /i)) — > Vq be a multiplicative map defined by L(X, fi) y-t 
L(A,/i)(A g) I n ). Its easy to see that M. is linear. First we show that M. is 
surjective. Let v ® Q(A,/i) be an arbitrary element of Vq. Construct L(A,/i) = 
AAi + + A3 where 

At = [v®A 20 v<g>An w®A 10 ] 
A 2 = [0 v®A 02 v®A m ] 
A 3 = [0 v <g> A 00 ] . 

Then obviously we have Ai EB A 2 EB A 3 = w <g> [A 20 A 02 A n A w A i A 00 ^ , 
so by Lemma l2~3l L(X. /i) is an A^-pre-image of v <g) Q(A, /i). The set of all M.- 
preimages of v <g) Q(A, /1) is L(A, /i) + KerAl, so all that remains is to compute 
KerAl. Further by Lemma 12.31 KerAl contains L{X,fi) = XA X + /j,A 2 + A 3 
that satisfies A x EB A 2 EB y4 3 = 0. The definition of "box- addition" implies that 
Ai,A 2 , A 3 are of the following form 

A = [0 -Y x -Z 1 ] 
A 2 = [Y t 0] 
A 3 = [Z x Z 2 0] 

where Y x , Z x , Z 2 G C 3nxn are arbitrary. This completes the proof. 



Example 2.6 In Example \2.4\ we achieve the linear two-parameter polynomial 
L(A,/u) G L(Q(A,/i)) by choosing 



-A 2Q 




-Am 




' " 


-A 00 + An 
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A 02 




-I + 2A 10 
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The standard linearization L(A, /x) G h(Q(X, fi)) is achieved by choosing 





~0 




'Aw 




~A 0l ~ 


Y 1 := 





,Z 1 := 





,Z 2 := 


-I 









-I 








Corollary 2.7 Dimension ofL(Q(X,fi)) = 9n 2 + 3. 

Remark 2.8 For quadratic one-parameter matrix polynomial 

P(A) = A 2 A 20 + AA 10 + A 00 , A i0 G C" x ™, i = 0, 1, 2, (12) 

a vector space Li(P) o/ matrix pencils of the form L(X) = X + AF G C 2rax2n is 
obtained in 121^ . Setting fi = m Q(X, fi) we have Q(X,0) = P(A). Then from 
the constructions of linear two-parameter polynomials given in Theorem \2.5\ it 
is easy to check that h(Q(X, fi)) = Li(P). In fact, if /j, — £/ien L(Q(A, //)) 
contains matrix pencils L(X) = XAi + A% G C 2nx2n where 

Ax = [v ® A 20 -Zi + v ® A 10 ] , A 3 = [Zi v <g A 00 ] , 

v G C 2 and Zi G C 2nxn zs arbitrary. Thus we obtain the same vector space 
of matrix pencils obtained in [21] for a given quadratic one-parameter matrix 
polynomial Q(X, 0) = X 2 A 20 + XA W + A 00 = P(A). 



2.1 Construction of linearizations 



It is not very clear that whether all linear two-parameter matrix polynomials 
in the space L(Q(A, / u)) are linearizations of Q(X,fi). For example, consider any 
L(X, /jl) G h(Q(X,/i)) corresponding to ansatz vector v — 0. Thus given a 
quadratic two-parameter matrix polynomial Q(X, fi) we need to identify which 
L(X, /jl) in h(Q(X, fj,)) are linearizations. 

We begin with a result concerning the special case of the ansatz vector 
v = ae\ where ei — [l 0] and ^ a G C. 

Theorem 2.9 Let Q(X, fi) = X 2 A 20 + XfiA u + /j 2 A 02 + XA 1Q + fiA 01 + A 00 be 
a quadratic two-parameter matrix polynomial with real or complex coefficient 
matrices of order n x n. Suppose L(X, fi) = XAi + fiA 2 + A 3 G h(Q(X, fi)) with 

ae% G 



3 , where 
An 



respect to the ansatz vector v 

Ax = [aei <g> A 20 —Y\ + ae\ 
A 2 = \Yi ae x <8> A 02 -Z 2 + ae x 
A 3 = [Zi Z 2 ae x ® A 00 ] , 



-Zi 
A i] 



A w ] 





~Yn 




~Z n 




Z\ 2 


Yi = 





,z x = 


Z 2 \ 


,z 2 = 


z 22 









Zzx 




z 32 _ 


L(X,fi) 


is a 


linearization 


ofQ{X,fj.) 





G C 3nxn , det 



Z 2 \ 
Z31 



z 22 

Z32 



^ 0. Then 
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Proof: By Theorem l2.5[ any linear two-parameter matrix polynomial L(X, fi) = 
XAi + fiA 2 + A3 E h(Q(X, fj,)) corresponding to the ansatz vector v = ae\ is of 
the form 



L(A,/i) 





aA 20 


-Y n + aA n 


-Zii + aA 


A 





-Y21 




—Z 2 \ 







-Y 31 




—Z31 




~Zn 


Z\2 OiA m 






+ 


Z 2 i 
Z31 


z 22 
Z32 
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Y u aA 02 
Y21 
I31 



-Z\i + aA 01 
—Z22 
—Z32 



Thus we have 

l(a,a«) 



fiY 21 + Z 21 
fj,Y 31 + Z31 



~XY 2 i + Z 22 
Z32 



W 3 (A,//) 
■A^2i — l^Z 22 

-XZ31 — flZ 32 



where W^A, fi) = aXA 20 + fiYn + Zu, W 2 (X, fx) = afiA 02 + XaAn 
and W3 (A, fi) = aXA w — XZn + afiA m — fiZ\ 2 + aA Q o. 
Define 

*I / 0" 

E{\fi) 



XY n + Z 



12 



a 

±1 7 

a 

ij 



Consequently, we have 



Q(X,fi) Wi(X,ii) 

^y 21 + z 21 
M y 31 + ^31 



W 2 (X,fi) 

"Ay^i + z 22 

^32 



Setting Y 21 = 



Y31 we have L(A, fi)E(X, fi) 

nx2n £ 



W(X,fi)= [Wi(X,fi) W 2 {\fi)] GC 
Since Z is nonsingular, we define 



Q(X,fi) W(X,fi) 
Z 

Z21 z 22 



where 



Z31 Z- 



32 



G C 



2nx2n 



Then we have 



F(X,li) 
F{X,Li)L(X,fj,)E(X,Li) 



I -W{X,fi)Z- 1 

z- 1 



Q(X,fx) 








Note that both E(X,fi) and F(X,fi) are unimodular polynomials. Hence we 
have detL(A,/i) = 7detQ(A, fi) for some nonzero 7 G C. Thus L(X, fi) is a 
linearization. This completes the proof. 

Let Q(X, fi) quadratic two-parameter matrix polynomial and L(A, fi) G L(Q(A, fi)) 
corresponding to an ansatz vector 7^ v G C 3 . Then the following is a procedure 
for determining a set of linearizations of Q(X, fi). 



Procedure to determine linearizations in L(Q(A,/i)): 
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1. Suppose Q(X,fx) is a quadratic two-parameter matrix polynomial and 
L(X, fx) = XAi + fxAi + A 3 G h(Q(X,/x)) corresponding to ansatz vec- 
tor v G C 3 i.e. L(X, fx) (A <g> 4) = u <g> Q(A, /i) . 



2. Select any nonsingular matrix M 



mn mi2 mi 3 
m 2 i m 22 m 23 
m 3 i m 32 m 33 



such that Mv 



ae\ G C 3 ,a 7^ 0. A list of nonsingular matrices M depending on the 
entries of v is given in the Appendix. 

3. Apply the corresponding block-transformation M <g> 4 to L(X,fx). Then 
we have L(X, fx) = (M £g> I n )L(X, fx) = XAi + fxA 2 + A 3 such that 



A x 

a7 2 

A 3 



cnei <g> A 20 — Fi + ae\ <g> A n —Z\ + ae\ ® A w 
Yi ae\ g) A 02 -Z 2 + aei <8> A i 
Zi Z 2 aei <g) A 00 



where 



















mnYu 








= (M§ 


5 4 


)Fi 


= (M$ 


5 4) 







m 2 iY 11 























_m 3 iYn_ 


















"Zh" 




muZu 


+ m 12 Z 21 


+ m 13 Z 31 




= {Md 


5 4 


)Z 1 


= (MS 




Z21 




m 2 iZ n 


+ m 22 Z 2 i 


+ m 23 Z 31 














Z 3 \_ 




_m 31 Z u 


+ m 32 Z 21 


+ m 33 Z 31 














Z\2 




m u Z 12 


+ mi 2 Z 22 


+ m 13 Z 32 




= (M§ 


5 4 


)Z 2 


= (M(S 


5 4) 


z 22 




m 21 Z 12 


+ m 22 Z 22 


+ m 23 Z 32 














z 32 




m 31 Z 12 


+ m 32 Z 22 


+ m 33 Z 32 



are arbitrary. 

4. For L(X, fx) to be linearization we need to choose Yi,Z±,Z 2 as follows. 
If 77i2i = 772 3 i = then choose Yu arbitrary; otherwise choose Y\\ = 0. 



for which 





~Z U 




Z\ 2 


Further we need to choose Z\ = 


Z 2 \ 


,z 2 = 


z 22 




z 31 




z 32 



det 



m 21 Z u + m 22 Z 2l + m 23 Z 3 i 
m 31 Z u + m 32 Z 21 + m 33 Z 3 i 



m 2 iZi2 + m 22 Z 22 + m 23 Z 32 
m 3l Z 12 + m 32 Z 22 + m 33 Z 32 



± 0. (13) 



From the construction of M given in the Appendix it is easy to check 
that we can always choose suitable Z\, Z 2 for which the condition (fl~3|) is 
satisfied. 
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3 Linearization of two-parameter quadratic eigen- 
value problem 

The quadratic two-parameter eigenvalue problem is concerned with finding a 
pair (A,|u)eCxC and nonzero vectors X{ G C Ui for which 

Qi(\,n) Xi = 0,i=l,2, (14) 

where 

Qi(\, //) = Ai + XBi + fid + A 2 A + XfiEi + fi 2 F t , (15) 

Ai, Bi, . . . , Fi G C niXni . The pair (A, //) is called an eigenvalue of flUJ) and x\ ® 
x 2 is called the corresponding eigenvector. The spectrum of a quadratic two- 
parameter eigenvalue problem is the set 

a Q ■= {(A, fi) G C x C : detQi(A, (i) = 0, i = 1, 2} . (16) 

In the generic case, we observe that (111]) has Ari\n 2 eigenvalues by using the 
following theorem. 

Theorem 3.1 (Bezout's Theorem, J^) Let f(x, y) = g(x, y) = be a system of 
two polynomial equations in two unknowns. If it has only finitely many common 
complex zeros (x, y) G CxC, then the number of those zeros is at most degree(f)- 
degree(g). 

The usual approach to solving ffT4"|) is to linearize it as a two-parameter 
eigenvalue problem given by 

/ L 1 {X, t i)w 1 = {A^) + XB( 1 ) + iJ l C^)w 1 = 

\ L 2 (\,fi)w 2 = (A^) + \B^ + fiC^)w 2 = {U) 

where A^,B^,C {i) G C m * xm %m t > 2m, i = 1,2, and w { = A ® x t . A pair 
(A,//) is called an eigenvalue of ( 1T7|) if Li(X,/J,)wi = for a nonzero vector u>j 
for z = 1,2, and W\ £g> 102 is the corresponding eigenvector. Thus the spectrum 
of the linearized two-parameter eigenvalue problem is given by 

a L ■= {(A, n) G C x C : detL^A, /jl) = 0, i = 1, 2} . (18) 

Therefore, in the generic case, the problem (JT7J) has m\m 2 > 477x772 eigenvalues. 

A standard approach to solve a two-parameter eigenvalue problem (|T7|) is 
by converting it into a coupled generalized eigenvalue problem given by 

Ai2 = XA z and A 2 z = /iA z (19) 

where z = w\® w 2 and 

A = ® C<® - C« ® £? (2 > 

A 2 = A« ® - B« <g> A( 2 ). 
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The two-parameter eigenvalue problem is called singular (resp. nonsingular) if 
A is singular (resp. nonsingular), see [21]. 

As mentioned earlier, we are interested in finding linear two-parameter poly- 
nomials Lj(A,/i) for a given quadratic two-parameter eigenvalue problem ( 1T41) 
such that gq = gl- Thus we have the following definition. 

Definition 3.2 Let ( fjjp be a quadratic two-parameter eigenvalue problem. A 
two-parameter eigenvalue problem [T7\ ) is said to be a linearization of [L~$ if 
Li(X, ii) is a linearization of Qi(X, Li). 

Thus if we consider a linearization of a quadratic two-parameter eigenvalue 
problem then gq = gl is guaranteed. It is also easy to observe that X\®Xi is an 
eigenvector corresponding to an eigenvalue (A, /x) of a quadratic two-parameter 
eigenvalue problem if and only if w\ ® W2 is an eigenvector corresponding to the 
eigenvalue (A, /i) of the linearization. 

Making use of the construction of linearizations for a two-parameter quadratic 
matrix polynomial described in section 2, we construct linearizations for a 
quadratic two-parameter eigenvalue problem. 

Theorem 3.3 Let [Lfy be a quadratic two-parameter eigenvalue problem. A 
class of linearizations of [Lfy is given by 



Li(A, Li)wi = (A (i) + XB® + LiC [i) )wi = 0, w t = A <g> x h i = 1, 2, 



where 



B {i) 
C (i) 



zf Zf a i e 1 ®A i 



n 



8> A 

CtiCi 



Ei 



-y/ + a iei 
Fi -Zf + a^i 



-Zf + a^i <g> Bt 

a 



0. 



(0 



rywi 








^12 





1 — 


^21 


j ^2 — 


^22 


_ _ 




. 31. 




7 (0 

. 32. 



7 (i) 

7 (i) 
^31 



^22 
7 (i) 
^32 



7^ 



Proof: Consider the linearizations Li (A, /i) = A^> + XB^' + fxC^ of Qj(A, //) 
associated with ansatz vector 7^ c^ei G C 3 ,z = 1,2 given by Theorem 12.91 
This completes the proof. 

Now we show that the linearizations for a quadratic two-parameter eigen- 
value problem described in Theorem 13.31 are singular linearizations. The follow- 
ing theorem plays an important role in the sequel. 



Theorem 3.4 The determinant of a block-triangular matrix is the product of 
the determinants of the diagonal blocks. 
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Now we have the following result. 



Theorem 3.5 The linearizations for (1$ ) derived in Theorem \3.3\ are singular 
linearizations. 



Proof: Consider the linearizations Lj(A, fi)wi 
0, i — 1, 2 of Qi(X, fi) where 



(A« + A5» + {j,C®)wi 



B 0) 








r(i) 



"^21 
"^31 



C (<) 



o 





OiiFi 





-^ 2 } + a % C t 

— ^22 
_Z, 32 



Consequently we have 

A = £? (1) ® C (2) - C« ® B< 2 > 



aiA g> C( 2 ) + ai^i) ® C (2) 






Y^®B^ a 1 F 1 ®BW 



Z 



(i) 

12 













{-Z$ + a 1 B 1 ) 

-z£> ® c& 

aid) ® 

£>(2) 

Zg } ® 5( 2 ) 



C(2) 



7(1) 
v 22 

Ki) 



Observe that Ao is a block-triangular matrix with one of the diagonal blocks is 
0. Hence by Theorem 13.41 we have detA = 0. This completes the proof. 

Remark 3.6 Note that given a quadratic two-parameter eigenvalue problem 
p^P we choose linearizations Lj(A, y) of Qi(\, n) associated with the ansatz 
vector § 7^ a^ei G C 3 ; and constructed linearizations Lj(A, fi)Wi = 0, Wi = A®Xj 
of However, we can derive a large class of singular linearizations by choos- 
ing linearizations Lj(A,/i) of Qi(\, /£) associated with ansatz vector 7^ G C 3 
described in section 2. 



4 Conclusions 

Given a quadratic two-parameter matrix polynomial Q(A,/i), we construct a 
vector space of linear two-parameter matrix polynomials and identify a set of 
linearizations of Q(\,fi). We also describe construction of each of these lin- 
earizations. Finally, using these linearizations we determine a class of singular 
linearizations for a quadratic two-parameter eigenvalue problem. 

Acknowledgement The author wishes to thank Bor Plestenjak and Rafikul 
Alam for their stimulating comments that have significantly improved the qual- 
ity of the results. Thanks are due to the referees for their constructive remarks. 
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Appendix 

Let O^aGC and e\ 



Given a vector v 



G C 3 we can always 



pick a nonsingular matrix M G C for which Mv = ae\ as follows. 



M = ( 



a/a 
l/a -l/b 
l/a -1/c 
b a/b ' 

-l/b 1/c 

1 
1 1 a/c 
1 1 

1 
'a/a 

1 
-l/a 1/c 

a/ a 0' 
1 
1 1 
a/a 1 
l/a -l/b 1 
-l/a l/b 

1 a/b 6 
10 
1 1 
a/a 
l/a -1/c 

1 



, if a^0,6^0,c^0 
if a = 0,6^0,c^0 
if a = 0,6 = 0,c ^ 
if a^0,6 = 0,c^0 
if a ^ 0,6 = 0,c = 
if a^0,6^0,c = 
if a = 0,6^0,c = 
if a^0,6 = 0,c^0. 
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